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ABSTRACT 

Recent XMM-Newton observations of the cooling flow in M87 indicate sharply decreasing oxygen, iron 
and silicon abundances within ~ 5 kpc of the galactic center. This result is unexpected since stellar mass 
loss and Type la supernovae are expected to produce pronounced central abundance maxima for all three 
elements. However, it has been suggested that many of the strong X-ray lines are optically thick and 
diffuse to larger radii in the cooling flow before escaping, falsifying the central abundances. We verify 
with radiation transfer calculations that this effect does indeed occur in the M87 cooling flow, but that 
it is insufficient to account for the M87 observations. We suggest that some source of continuous opacity 
is required to reduce the central X-ray line emission, perhaps by warm gas at T ~ 10 5 — 10 6 K. The 
radial surface brightness profiles of X-ray resonance lines are also sensitive to turbulence in cooling flows 
which reduces the line center optical depths considerably. Turbulence may provide sufficient energy to 
continuously heat the warm absorbing gas. 

Subject headings: galaxies: elliptical and lenticular - galaxies: cooling flows - galaxies: interstellar 
medium - X-rays: galaxies - X-rays: galaxy clusters 



1. DIFFUSION OF X-RAY LINE PHOTONS IN M87 

This estimate of the diffusion of X-ray line photons in 
M87 has been motivated by the exceptionally low abun- 
dances of iron, silicon and oxygen within ~ 5 kpc of the 
center of M87 as recently observed by Bohringer et al 
(2000) with XMM-Newton. We explore the possibility 
that these abundance variations are artifacts of the out- 
ward diffusion of optically thick X-ray lines as suggested by 
Gil'franov et al. (1987), Tawara et al. (1997), Shigeyama 
(1998) and Bohringer et al. (2000). We find that the ex- 
pected outward migration of line photons does occur, but 
is insufficient to account for the central photon depletion 
observed. Instead, an additional source of continuous ab- 
sorption is indicated, as suggested by Buote (2000), and 
we speculate about its physical nature. Before discussing 
our results, we begin with a review of the properties of the 
M87 cooling flow and describe our line transfer calculation. 

2. VARIATION OF GAS TEMPERATURE AND DENSITY 

The electron density distribution in the M87 cooling flow 
is taken from Nulsen and Bohringer (1995), corrected to 
distance 17 Mpc (1' = 4.945 kpc): 

n e , g = n [l + {r/r o )r>}- 1 (1) 

with no ~ 0.120 cm~ 3 , p — 1.18 and r Q = 0.5r e where 
r e = 8.625 kpc is the effective radius. According to 
Bohringer et al. (2000) these observations underestimate 
the central density which should instead be represented by 
a (3 model fit : 

n e , /3 = n 1 [l + (r/r /3 ) 2 ]- 3/3/2 (2) 
4 Chandra Fellow. 



with ni = 0.35 cm' 3 , (3 = 0.47 and rp = 1.97 kpc. This 
inner density distribution is based on ROSAT HRI obser- 
vations reported in Bohringer (1999). Since the central 
X-ray emission in M87 is known to be partly non-thermal 
(and time variable) , this latter density determination may 
be uncertain. Nevertheless, in the following we assume 
that n e {r) is given by the maximum of n ejS (r) and n e ^{r) 
which intersect at 11.9 kpc. 

The (single-phase) gas temperature distribution is found 
by combining the results of Bohringer et al. (2000) 
within 65 kpc with more distant temperatures measured 
by Nulsen & Bohringer (1995). These results can be fit 
with the fitting function: 

T = To( r -) L-T* +(-!)*] 1 (3) 

\r tj l(r + r ot ) \r m ) 

withT = 1.5xl0 7 K, r m = 5.695, r ot = 1.53 and q = 0.22; 
all radii are in units of r e . This temperature profile rises 
from the center toward a broad maximum, similar to ther- 
mal profiles in other well-observed bright elliptical galaxy 
cooling flows (Brighenti & Mathews 1997). The location 
of the temperature maximum in M87 is at an unusually 
large radius presumably due to the relatively high ambient 
temperature (2-3 keV) in the Virgo cluster. The central 
cooling of the gas can be understood as the mixing of hot 
inflowing gas from the Virgo cluster with gas lost from 
evolving stars in M87 which orbit with a lower virial tem- 
perature than that of the surrounding cluster (Brighenti 
& Mathews 1999a). The characteristic deep central tem- 
perature minimum in M87 indicates that most of the gas 
in M87 has not been substantially heated by the jet and 
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radio source. 

3. ABUNDANCES, LINE EMISSION AND ABSORPTION 

We assume meteoritic solar abundances with log(a,g) = 
log(ns/nff) = —3.07 and —4.45 for O and Si respectively. 
The radial abundance variation of each clement in the hot 
gas relative to solar meteoritic is represented with A E (r). 
The line emissivity is assumed to vary with local temper- 
ature according to the emission coefficients determined 
by Landini & Monsignori Fossi (L-MF; 1990): 



Si = A E (r)n e n H aAe(T) (erg cm 



3 s- 1 ) 



where a is a factor of order unity that corrects L-MF abun- 
dances to solar meteoritic and nu — n e {& — 3fi)/(2 + //) is 
the local proton density in a plasma with molecular weight 
/i = 0.61. 

The emission lines that we consider arise from the 
ground state and are Doppler broadened by thermal and 
turbulent velocities. The absorption coefficient is 

k(x) = n H Xia E A E ( ^—f gu — eT^ (cm -1 ). 
\m e c 7T ' Av D ) 

Here a E A E (r) is the local abundance of the line-emitting 
element and x% (r) is the fractional abundance of the line- 
emitting ion interpolated from the results of Sutherland & 
Dopita (1993) for collisional ionization equilibrium. Pho- 
tons are assumed to be emitted in a Gaussian profile cen- 
tered on the line energy E Q — hv and 
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is the dimcnsionlcss line frequency with 
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where A Ewt is the atomic weight of element E. The line 
center optical depth in the radial direction is found by 
integrating dr = n{Q)dr. 

We are taking Xi(T) and Ae(T) from two different 
sources. However, the error from this possible inconsis- 
tency is unlikely to be large since collisional ionization 
fractions Xi(T) found by most authors, including L-MF, 
are in good agreement with those of Arnaud & Rothenflug 
(1985). 

4. MONTE CARLO PROCEDURE 

To estimate the diffusion of X-ray line photons, we rep- 
resent the inner cooling flow of M87 with N = 50 radial 
zones, each having different but individually uniform tem- 
perature and density. The zone spacing r n is chosen so 
that the line luminosity is the same from each zone. The 
Monte Carlo process begins by choosing the zone of origin 
of each photon based on equal weighting among the zones. 
The exact birth radius of each photon within the chosen 
zone is then found from 

r b = [rl+Tl(rl +1 -rl)]^ 

where 1Z is a random number between and 1. The ini- 
tial direction of the photon is w — cos 9 — 21Z — 1 where 
1Z is another random number and 6 is the angle between 
the photon velocity and the direction to the center of the 
cooling flow. The initial frequency x is chosen by random 



selection from a Gaussian probability density. The dis- 
tance that the photon moves from its origin r a to the local 
zone boundary r z is 

£ = wr + a[r 2 z -r 2 {l-w 2 )}^ 2 . (5) 
Here a = 1 and r z = r n+ \ for outward moving photons 
(w < w cr ) and a = — 1 and r z = r n for inward moving 
photons (w > w cr ) where w cr = [1 — (fTi/fo) 2 ] 1 / 2 . In the 
zone of photon birth r Q — rj,. If the photon continues to 
move through several zones without scattering, w and r a 
remain constant and Equation (5) can be used to find the 
total path length moved L. The distance £ traversed by 
the photon across each zone must be compared to the free 
path in the zone: 

l s = - \n{Tl)/K n {x n ) 
where K n (x n ) is the absorption coefficient of the pho- 
ton in zone n. Since Avb varies from zone to zone it 
is necessary to rescale the frequency to the local zone n 
(r„ < r < r„+i): x n (r) = x n {r b )[Au D (r b )/Au D (r)}^ 2 . If 
i s > £, the photon moves into the next zone, but if £ s < I 
the photon is scattered at radius 

^r s = {r 2 + Ll-2r L s w) 1 l 2 
where L s = ^2£ n + £ s - The photon is then scattered into 
a new random direction beginning a new trajectory with 
r a = r s , assuming complete redistribution in the Gaus- 
sian line profile. The scattering is conservative so all pho- 
tons must ultimately flow out beyond the outermost zone. 
When the photon reaches the outer boundary of the gas 
ft = riV) where the optical depths are very small, it is 
assumed to escape freely. 

5. RESULTS AND DISCUSSION 

We describe transfer calculations for two lines: SiXIV 
2.003 keV (La) and OVIII 0.652 keV (La), both having re- 
values f gu = 0.420. Since these resonance lines are among 
the most optically thick lines of Si and O at M87 gas tem- 
peratures, the apparent abundance distribution inferred 
from each of these lines - subsequent to spatial diffusion - 
must correspond to the largest possible central abundance 
depiction. For each line-emitting element the apparent ra- 
dial abundance profiles observed with XMM-Newton can 
be fit with 

qob 

/ 

+ 
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Tmob 



Totob 



f*mob 



(r + r ot ob) \r mob 



( 6 ) 

The apparent (meteoritic) abundance dependence of Si 
and O are approximated with: (A ob , r mo b, r oto b, qob) = 
(0.1 , 7.4, 0.34, 0.7) and (0.05, 8.0, 0.5, 0.94) respectively; 
Tmob and r oto b are in kpc. The line center optical depths 
in the three lines evaluated with A^ b \r) (Fig. 1) are 
smaller than those suggested by Bohringer et al (2000) 
in part because we use the higher temperatures observed 
instead of their isothermal approximation. If the observed 
abundance maxima are due to outward photon diffusion 
by conservative scattering before escape, the total X-ray 
line luminosity within some larger radius (where r <C 1) is 
unaffected by the radiative transfer. 

According to this line diffusion hypothesis, the true ra- 
dial abundance profiles are expected to decrease monoton- 
ically as we represent with 
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The first term accounts for the local enrichment by stel- 
lar mass loss and SNIa explosions, and the second term 
represents the abundance of cluster cooling flow gas flow- 
ing into M87, but the parameterized fit is not intended 
to exactly quantify these two enrichment sources. The 
coefficients in Equation (7) are chosen so that the total 
line luminosities within 50 kpc are very close to luminosi- 
ties based on A^ (see below): (A Q , rs(kpc), c, A c f) = 
(1.0, 35, 1.4, 0.15) and (0.4, 12, 2, 0.14) for Si and O re- 
spectively. The A c f are chosen for agreement with the 
XMM-Newton abundances at r > 30 kpc; te is assumed 
to be similar to r e = 8.625 kpc appropriate for a stellar 
contribution; c is chosen for a reasonable fit to the XMM- 
Newton data in 10 < r < 30 kpc; A Q is adjusted until 
the total luminosity within r = 50 kpc (where all lines are 
optically thin) matches that of each line computed with 

Ag b \r). With these assumed intrinsic abundance profiles 
Ae(t), the total line luminosities within 50 kpc exceed the 
approximate observed luminosities of the Si and O lines by 
about 2 and 13 percent respectively. We exceed the oxy- 
gen line luminosity based on Aq^ (r) since this abundance 
is less accurately determined by the XMM-Newton data. 
The cumulative line center optical depths are somewhat 
greater when evaluated with Ae(t) than with Ag b \r) as 
shown in Figure 1. 

Monte Carlo radiative transfer for each line has been 
calculated with 10 4 photons, 50 spatial zones and abun- 
dances Ae{t). Results are illustrated in Figure 2. The 
light dotted histograms show the abundance of each ele- 
ment that would be inferred by observing each line in the 
optically thin limit, Asibirth) = AE(r)N p h,bl ' {N) , where 
N p h,b is the number of photons created in each zone and 
(N) = 200 is the mean number of photons per zone. The 
adjacent light-line histograms show the variation of the 
apparent abundance indicated by each line after radiative 
diffusion, A^(escape) = AE(r)N p h, e / (N) , where N p h if . is 
the number of photons that escape from each zone with no 
further scattering. Ae (escape) can be compared directly 
with the observed abundances. 

The light solid-line histograms clearly indicate central 
photon depletions or flattening within ~ 5 kpc. The ap- 
parent Si abundance based on the 2.003 keV line is clearly 
self-absorbed and matches the observations quite well at 
r > 3 kpc. But the amount of photon depletion with con- 
servative scattering is clearly insufficient to explain the 
observed O abundance. It is not possible to deepen the 
central minima by increasing the O abundance Ao (r) (and 
line opacity) in r < 5 kpc since the global line luminosities 
would exceed those observed. The OVIII line luminos- 
ity may already be too high by about 13 percent. The 
central oxygen line luminosity (abundance) only levels off 
and does not decrease toward the origin as observed. Ev- 
idently line photons are being lost by absorption. The 
central abundances of Si and O - Ae(0) — 1.15 and 0.54 
respectively - cannot be trusted if the line scattering is 
non-conservative. Since the continuous opacity in a warm 
plasma with T w < 10 6 K decreases rapidly with photon en- 
ergy, a oc E~ 3 , the OVIII line should suffer ~ 30 times 
more absorption than the SiXIV line, in qualitative agree- 
ment with the profiles in Figure 2. 

Shigeyama (1998) noted that if the outward diffusion of 



strong X-ray lines is not considered in the data reduction, 
the central gas density can be underestimated by ~ 20 — 50 
percent. To explore this, we repeated the calculations with 
n £j/ 3 increased by 1.5. Our conclusions are strengthened 
by this adjustment in the gas density. With higher n e ^ 
the luminosity of the lines within 50 kpc can be matched 
with a slightly lower central abundance, resulting in cen- 
tral minima in Ae (escape) (r) even less pronounced than 
those in Figure 2. 

Studies of optical line emission in the central regions of 
galactic and galaxy group cooling flows reveal widespread 
random motions typically about 0.2 — 0.4 of the sound 
speed in the hot gas (Heckman et al. 1989; Caon et al. 
2000). Since small, optically visible line-emitting regions 
at T ~ 10 4 K are likely to be strongly coupled to the ambi- 
ent hot gas, we infer that the hot gas must share the same 
turbulent velocities. Turbulence is expected to strongly 
influence the spatial distribution of X-ray resonance lines, 
particularly for more massive line-emitting ions (Equation 
4) . To demonstrate this, we repeated the transfer calcula- 
tions for both lines with v^ urh — A4 2 ui , & (5^4,e/6/z) assuming 
a turbulent Mach number M tur b = 0.3. In the presence 
of turbulence, which we regard as more realistic, the line 
center optical depths are greatly reduced (Fig. 1), there is 
less outward diffusion (Fig. 2), and the apparent central 
depletion of silicon photons is significantly reduced. 

We have also computed line profiles for the FeXXV 6.69 
keV (He 4) resonance line discussed by Bohringer et al. 
but have not discussed its profile here because emission 
in this line is unlikely to be representative of the intrinsic 
iron abundance. The FeXXV 6.69 keV line is very sensi- 
tive to the (low) cooling flow temperatures T(r) in r < 10 
kpc where the line transfer is most relevant (e.g. Fig. 
2 of Buote, Canizares & Fabian 1999) and relatively few 
FeXXV line photons are created in this low temperature 
region. Because of this, the x 2 iron abundance determi- 
nation in M87 by Bohringer et al. is dominated by the 
multitude of strong Fe L lines near ~ 1 keV. The many 
lines that constitute the Fe L feature have different atomic 
parameters, so no single L line can faithfully represent 
the iron abundance, as we have assumed for the O and 
Si Lyman-a lines. Since the continuous plasma opacity 
varies as a oc E~ 3 , the opacity at the FeXXV 6.69 keV 
is ~ 1000 times less than the opacity at the OVIII 0.652 
keV line. An observation of the equivalent width of the es- 
sentially unabsorbed FeXXV 6.69 keV line with projected 
radius in M87 would allow a calculation of the intrinsic 
iron abundance distribution, provided T(r) is sufficiently 
well determined. 

We conclude that the central minima in apparent abun- 
dances and line luminosities in M87 cannot be understood 
solely in terms of conservative photon scattering. An ad- 
ditional source of continuous opacity seems to be required. 
We speculate that the source of this opacity is the warm 
absorbing gas at T - 10 5 - 10 6 K as suggested for M87 
by Buote (2000). Buote's conclusion is derived from mod- 
eling the ROSAT PSPC data of M87 with coronal plasma 
emission modified by a single absorption edge and Galactic 
absorption. The rest energy of the edge is determined to lie 
at ~ 0.6 keV, near the position of the OVIII La line. The 
absorption derived from the ROSAT data (r ~ 0.5 — 1) is 
spatially more extended than that required to account for 
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central minima in the X-ray resonance lines. Because of 
the limited energy resolution of ROSAT, the effect of using 
a simple absorption edge cannot be distinguished from a 
model in which the O abundance is allowed to take low 
values as done by Bohringer et al (2000). 

The central declines in Fe and Si abundances in the 
Bohringer ct al. data are difficult to understand since we 
expect hot gas closer to the galactic core to be more en- 
riched by centrally peaked stellar mass loss and Type la 
supernovae (Brighcnti & Mathews 1999b). Since the oxy- 
gen abundance found by Bohringer et al. has a generally 
negative gradient beyond ~ r e , it is apparent that the 
stellar and Type la ejecta are more oxygen-rich than cool- 
ing flow gas entering M87 from larger radii, so a central 
peak in the O abundance would also be expected. The 
need for an additional source of continuous opacity in our 
line transfer calculations lends new support to the absorp- 
tion interpretation. If this absorption arises from warm 
gas, then its mass is similar to that of the hot gas (Buote 
2000). 

The energy source that heats the warm gas is of great 
interest. While some of the cooling flow gas may be heated 
by the M87 jet, this region is avoided in the Bohringer et 
al. observations; the deep central minimum in their T(r) 
shows no evidence of AGN heating. The most likely energy 
source for heating the warm gas is the turbulent motion 
of the hot gas as indicated by the dispersion velocity of 



optical emission lines Vt ~ 200 km s _1 in M87 (Heckman 
et al. 1989). These velocities are subsonic in the hot gas, 
but would be supersonic in the warm, X-ray absorbing 
gas. Colliding regions of warm gas would be shock-heated 
at at a rate H = fn c m p vf jr c ergs cm~ 3 s _1 where / is 
the filling factor of warm gas clouds of mean radius r c and 
electron density n c . Although we expect that some net 
cooling occurs, if shock heating balances radiative cooling 
C = n 2 A ergs cm' 3 s' 1 then the cloud size can be esti- 
mated: r c ~ fm p vf/n c A ~ 400/ pc for v t — 200 km s _1 , 
n c = 0.1 cm -3 and cooling coefficient A = 10~ 22 erg cm 3 
s _1 . With / ~ 0.5, this cloud size is consistent with the 
scale of irregularities observed in the optical line emission 
and the relative uniformity of the warm plasma absorp- 
tion. 

Finally, we speculate further that such extended absorb- 
ing regions are generally present in cooling flows and ac- 
count for the apparent deficiency of emission lines charac- 
teristic of cooling ions in the centers of these flows (e.g. 
Peterson et al. 2000; McNamara et al. 2000). 
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Fig. 1. — Variation of line center optical depths in the central cooling flow of M87 in two lines: SiXIV 2.003 keV and OVIII 0.652 keV. 
The solid and dashed lines show optical depths resulting from the model Ae(t) and directly from the observations A^ b \r) respectively. The 
heavy dashed lines are line center optical depths based on Ae(t) but in the presence of turbulence with M tur i, = 0.3. 

Fig. 2. — Radial variation of silicon and oxygen abundances in the M87 cooling flow. The observed points including symbol type have 
been taken directly from Bohringer ct al. (2000). The smooth lines passing through the observations are A^ b \r) found from Equation (6). 
The dotted histograms show the variation of the elemental abundances Ag (birth) based on each line using the Monte Carlo program in the 
optically thin limit; apart from statistical fluctuations this histogram is identical to A^fr). The light solid histograms are the abundance 
profiles Ae (escape) that would be inferred from each resonance line after Monte Carlo spatial diffusion; the line center optical depths for 
these calculations are the light solid lines from Figure 1. The heavy histograms show the effect of turbulence with M tur f, = 0.3. 
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